Mott transition in the asymmetric Hubbard model at half-filling within dynamical 

mean-field theory 
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We apply the approximate analytic methods to the investigation of the band structure of the 
asymmetric Hubbard model where the chemical potentials and electron transfer parameters depend 
on the electron spin (type of quasiparticles) . The Hubbard-I and alloy-analogy approximations are 
the simplest approximations which are used. Within the alloy-analogy approximation, the energy 
band of particles does not depend on the transfer parameter of particles of another sort. It means 
that the gap in the spectrum opens at the critical value U c that is the same in two different limiting 
cases: the Falicov-Kimball model and the standard Hubbard model. The approximate analytic 
scheme of the dynamical mean-field theory is developed to include into the theory the scattering 
of particles responsible for the additional mechanism (due to the transfer of particles of another 
sort) of the band formation. We use the so-called GH3 approach that is a generalization of the 
Hubbard-III approximation. The approach describes the continuous Mott transition with the U c 
value dependent on a ratio of transfer parameters of different particles. 

PACS numbers: 71.10.Fd, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

The asymmetric Hubbard model is considered as a gen- 
eralization of the Falicov-Kimball model Q and the Hub- 
bard model Q ■ This model describes a system with two 
sorts of mobile particles (ions, electrons or quasiparti- 
cles) with different hopping integrals and different values 
of chemical potentials. 

The hopping of particles is described by creation and 
annihilation operators and transfer parameters tfj. The 
Hamiltonian of the model is 



H 



H, 



ija 



where the singe-site part 



Urii 



nil 



(1) 



(2) 



includes chemical potentials n a and a local on-site repul- 
sion U {n ia = a\ a a ia ). 

Due to a complexity of the problem the asymmetric 
Hubbard model is much less investigated than its lim- 
iting cases: the Falicov-Kimball and Hubbard models. 
The model was investigated in the large U limit at half- 
filling. The effective Hamiltonian corresponding to the 
anisotropic Heisenberg model was derived H 3 and the 
effective antiferromagnetic interaction was discussed 
A phase separation phenomenon was considered in the 
ground state of the asymmetric Hubbard model @ . Some 
progress_has_been achieved in the case of one-dimensional 
chains 
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The essential achievement in the theory of the strongly 
correlated electron systems is connected with the devel- 
opment of the dynamical mean-field theory (DMFT) pro- 
posed by Metzner and Vollhardt ^3 f° r the Hubbard 
model (see also Ref. |ll| for a review). This method is 
based on the assumption that the self-energy is a pure 
local (a single-site) quantity (w) = Y, a {u))5ij and it be- 
comes exact in the limit of infinite dimensions 12]. The 
current growing interest to DMFT is related to the new 
developing technique combining DMFT and the density 
functional theory within the local density app roximation 
(LDA). The new approach LDA+DMFT [lj| allows one 
to calculate the electronic structure of real materials cor- 
rectly taking into account the strong local correlations 
(see Refs. 0, for reviews). 

The Falicov-Kimball and Hubbard models were inten- 
sively investigated within DMFT (see reviews Refs. 0, 
16). In the Falicov-Kimball mo del, par ticle densities of 
states can be calculated exactly [13 111 El El El El . 
However, the spectrum of localized particles was in- 
vestigated mostly at half-filling. In the case of the 
half-filled Hubbard model, a number of numerical and 
analytic approximate methods was developed and ap- 
plied to the investigation of metal-insulator transitions 

iMElElElElilElElEiElEi. 

We use the approximate analytic approach to the in- 
vestigation of the band structure of the asymmetric Hub- 
bard model in DMFT l34l . This method was developed 
for the Hubbard model |35ll36| and it is based on a map- 
ping of the problem onto the effective single-site Hamil- 
tonian with the auxiliary Fermi-operators describing the 
environment of a given site. The approach is based on the 
equations of motion for Hubbard operators followed by 
the different time decoupling of the higher order Green's 
functions. The irreducible parts are separated off us- 
ing projecting on the basis of fermionic Hubbard oper- 



2 



ators. This approach gives DMFT equations in the ap- 
proximation that is a generalization of the Hubbard-III 
approximation and includes as simple specific cases the 
modified alloy-analogy approximation |37j and the Hub- 
bard-III approximation |38f. We call it the generalized 
Hubbard-III (GH3) approximation. Recently, an alter- 
native decoupling scheme in an equation of motion ap- 
proach to the Hubbard model in infinite dimensions has 
been proposed in Ref. 

This paper is organized as follows. In Section [nj we 
review the formalism of DMFT, where particle hopping 
can be introduced in two different ways: using the coher- 
ent potential or the effective single-site Hamiltonian with 
the auxiliary Fermi-filed . The equation of motion ap- 
proach with the projecting technique and the different- 
time decoupling scheme is described in Section IlIII This 
approach gives the DMFT equations within the gener- 
alized Hubbard-III approximation. In Section IIVI sim- 
pler approximations to the asymmetric Hubbard model 
are introduced. In Section [V] some analytic properties 
concerning continuous metal-insulator transitions are de- 
rived. Our results are discussed in Section IVII and con- 
cluding remarks are given in Section IVIII 



ffmt(r) = f At' J2 Mr r')4 (r)a CT (r'), (10) 



where 



Ho = H i: 



(11) 



and J ct (t— t') is the coherent potential. This function de- 
scribes the propagation of a particle in the environment 
without going through the given site between moments r 
and t'. It is determined self-consistently from the condi- 
tion that the same irreducible part determines the lattice 
single-site Green's function 

G„(u, n ) = G^K) = ^]TG CT K,k), (12) 



G a (uj n ,k) 



(13) 



and the Green's function of the effective single-site prob- 
lem 



II. DYNAMICAL MEAN FIELD THEORY 

The particle Green's function in the Matsubara repre- 
sentation is defined as follows 



(3) 



S CT (uj n ) - Jcr(Wn) 



(14) 



To investigate the Green's functions in the time rep- 
resentation the analytic continuation from the imaginary 
to real axis (iuj n — > u> + ie) is performed: 



G CT K) - G a {u>) = 2ir{{a a \at)) w 



(15) 



G£K) = / G£(r)e- 
Jo 



r dr, 



w(2n+ 1) 





(4) 



where T T denotes imaginary time ordering, and (5 = 1/T 
is the inverse temperature. 

For the Hamiltonian £[]), the particle Green's function 
can be written as a solution of the Larkin equation |40|: 

GijM = Sij CT (w„) + ^KlCCjK), (5) 



or in momentum space 

Go.(a; n ,k) = E a (u> n ,k) + 3 CT (w„, k)i£G CT (w n , k), (6) 

where 2 (T (a; Tl ,k) is the total irreducible part; it is con- 
nected to the Dyson self-energy by the relation: 

S^ 1 ^, k) = iuj n + n a - T, a (uj n , k). (7) 

In the limit of high lattice dimensions d — > oo, the 
irreducible part becomes a single-site quantity |l2t l4l| : 



^ija(,^n) — ^a^ij: ^cr(^'n;k) — ^ a (uj ri 



(8) 



The function S,j(a;„) or S ff (a; n ) is calculated using the 
auxiliary single-site problem. This problem corresponds 
to the following replacement 



-0H 



- / drff int (r) 
Jo 



(9) 



The self-consistency condition is rewritten in the fol- 
lowing form 



Pl(t)dt 



G a {w) 



(16) 
(17) 



where (|13|) is inserted into (|12fl and the sum over k is 
replaced by the integration with the noninteracting den- 
sity of states. The infinite dimensional hypercubic lattice 
and the Bethe lattice are the most popular lattices inves- 
tigated in the dynamical mean-field theory. A Gaussian 
noninteracting density of states corresponds to the hy- 
percubic lattice with nearest-neighbour hopping: 



t a \pK 



exp 



£ 



(18) 



where the transfer parameters must be scaled to keep 
the kinetic and interaction energies of the same order of 
magnitude in the infinite dimensional limit [Toj ]: 



(19) 



2Vd 



when i and j are nearest-neighbours and zero otherwise. 
In the case of the Bethe lattice with the coordination 



3 



number z — > oo, the noninteracting density of states is 
semielliptic 

1 



P T hc (e) = ^VK^, \e\<2U, (20) 

and the hopping integral for nearest-neighbours is scaled 
as follows: 



t = 



(21) 



and the two-time Green's function G CT (w) = 27r((a .|aj T )) w 
is expressed as: 

G a = 2n[((X Q °\X°°)) UJ +(((X a °\X 2 *)) UJ 

+ C{{X S2 \X M )) U + {{X S2 \X 2S )) U \ , (27) 

where the following notations for sort indices are used: 
(j = B, ( = + for a = A and a = A, £ = — for a = B. 

To calculate the Green's functions in (|27|) . we use the 
equations of motion for Hubbard operators: 



Besides the hypercubic and Bethe lattices, the nonsym- 
metric densities of states can be considered. The infinite 
dimensional generalization of a face-centred-cubic lattice 
that was used to describe ferromagnetic order in the Hub- 
bard model 01 i s an example of such a lattice. 

A solution of the single-site problem © gives the par- 
ticle Green's function G a (ui) as a functional of the coher- 
ent potential J a (uj). This dependence supplemented by 
the self-consistency condition forms a closed set of equa- 
tions for the single-site Green's function and the coherent 
potential. 

We reformulate the single-site problem introducing the 
effective Hamiltonian as it was done in Ref. |35f 



(22) 



The auxiliary fermionic operators £ CT are introduced 
to describe particle hopping between the site (H ) and ef- 
fective environment defined by the Hamiltonian H$ . The 
coherent potential is expressed as the Green's function in 
unperturbed Hamiltonian Hg: 



(23) 



Thus, an explicit form of H% is not required to solve the 
problem. 



dt 



(28) 



The commutators (12811 are projected on the subspace 
formed by fermionic operators X 0<7 and X" 2 : 



[X\H cS \ =aJX°° 



cP 2 X* 2 



Z* 1 



(29) 



The operators Z "^ 2 ^ are defined as orthogonal to the 
operators from the basic subspace |34L l35l l3q : 



(30) 



/fz 0a i a2 ) ; x "^}) = 0. 



Thus, these equations (|30|l determine the projecting co- 
efficients a^^ 2 ^ which are expressed in terms of mean 
values (£ CT X P9 ) combined into the following constant 



^ = (^x*°) + ((x° 2 t;i) 

and mean values of the Hubbard operators 



A pq = {X" + X™), A 0a = l- 



A 2S = 



(31) 



(32) 



This procedure leads to the Green's functions of the 
^o<T(a2)|j£ CT o(2cr)^ ty p e . The similar procedure can be 
applied with respect to the second time argument. As a 
result, we come to the relations between the components 
of the Green's function G a and scattering matrix P a . In 
a matrix representation, we have 



III. GENERALIZED HUBBARD-III 
APPROXIMATION 

We use the single-site Hamiltonian of the asymmetric 
Hubbard model written in terms of Hubbard operators 

H Q = - [^(X™ + X 22 )] + UX 22 , (24) 



acting on the basis of single-site states \ua, tib) 
|0) = |0,0), |A) = |1,0), 



|2) = |1,1) 



B = 0,1 



(25) 



In this case, the particle creation and annihilation oper- 
ators are expressed as 



a a =X 0a + CX S2 , 



(26) 



where 



2tt 



G a — Gq + Gq P a Go , 



( ((X 0a \X° )) {{X 0a \X 2s )) 



y {{x s2 \x M )) {{x s2 \x 2s )) 

and nonperturbed Green's function Gq is 



oj — b. 



'A 0rJ 
A 2S 



where 



D a = (uj - a a )(uj - b a ) - 



V 2 



A aA 2 , 



Va 



(33) 



(34) 



(35) 



(36) 



U-n a + -f-<p a . (37) 
A 2S 
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The scattering matrix 



Pa = 27T 



A 



Off 







A 2S 



((Z°°\Z° )) {{Z 0a \Z 2s )) 
{{Z 92 \Z a0 )) {{Z a2 \Z 2s )) 



A 



Oct u 

A- 



(38) 



contains the scattering corrections of the second and the 
higher orders in powers of V a . The separation of the 
irreducible parts in P a enables us to obtain the mass 
operator M a = P a \ii and the single-site Green's function 
expressed as a solution of the Dyson equation 



(1 - gim„)- x gi 



(39) 



The mass operator M„ is calculated in zero approxi- 
mation 



M„ = P (0) 



(40) 



where the time correlation functions related by the spec- 
tral theorem to the irreducible Green's functions are cal- 
culated using procedure of different-time decoupling. In 
our case it means an independent averaging of the prod- 
ucts of X and £ operators. For example, 

(C (t)(x m + x™)t(x oa + x-)O ir 

« ((X m + Xnt(X QQ + X™))(£(t)U- (41) 

For simplicity we calculate the correlation functions for 
Hubbard operators in zero approximation 



(pc uu + X aa ) t {X m + X aa )) 



■X c 



= A 



Oct- 



(42) 



In this example, the following Green's function is recon- 
structed from the correlation functions using the spectral 
representation pi I3I IH 143] : 



I(u;)=AoA(UC))i = 



A 



Oct 



2irV 2 



Ja(L0). 



(43) 



Using the above procedure, we can obtain the final 
expressions for the total irreducible part: 



A 



Oct 



A 2 a 



u + [i a — fi CT (w) uj + \i a - u — ricr(w) 

t-fi„(w), 



(44) 



where 



and 



Q a (u>) = J a {w) 



R*{oj) 



(X c 



x c 



-J S (UJ + \l a - fig) 



(45) 



{X aa - X s 
2^ 



dw'Im J s (u)' + i0+) (3uj' 
uj - uj' — \i s + \x a 2 



(X 00 + X 22 ) . 

H 7. Ja{U - [M a - Mct - w ) 



-<*°°-* 22 > /■d*W,(-a/-iO+) tanh /^ 

27T J u> - u>' + fl & + [i a - U 2 

—00 

The average values of the Hubbard operators are de- 
termined using the spectral representation of correspond- 
ing Green's functions. The particle density of states is 
calculated as an imaginary part of the interacting single- 
particle Green's function: 

p a (uj) = -- lim lmG a {uj + it]) (47) 

7T r;^0+ 



giving the concentration 



-f 00 



p a (m)duj 
e f3uj + 1 



(48) 



The Green's functions ((t; a \X pq }} u that are required to 
find ip a can be calculated from the exact relation 

v a ((Ux pq ))u> = JA")(Mx pq )h (49) 

derived in Ref. Thus, the approximation given by the 
equations l|44 [l -l|46 )l we called the GH3 approximation. 

We consider below the special case at half-filling (fj,A — 
Hb = U/2, riA = Ub = 1/2), and due to the particle-hole 
symmetry we have: 



<Ax = 0, 



Q a (u) = J a {uj) + 2J CT (w), 



and 



U 2 



4[uj-J a (oj)-2J„(uj)] 



(50) 
(51) 



(52) 



For the standard Hubbard model (pA — is, J a — Jb) 
the approximate solution of the single-site problem gives 
the usual Hubbard-III approximation. 

In the case of the Falicov-Kimball model (t B = 0, 
Jb{<uj) = 0) at half-filling, GH3 gives the following ex- 
pression: 



G b {uj) 



uj — 2Ja(uj) 



uj 2 - [/ 2 /4-2wJ A (w)' 



(53) 



where the coherent potential of itinerant particles J a cor- 
responding to the result of the alloy-analogy approxima- 
tion is calculated exactly. 



IV. OTHER APPROXIMATIONS TO THE 
ASYMMETRIC HUBBARD MODEL 

A. Hubbard-I approximation 

It was the Hubbard-I approximation 2] which in the 
simplest way described band forming in the Hubbard 
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Pa 



0.2 



0.0 



(C)n B =0.2 


AA 




H-l 







(D)n B =0.2 
U=3_ 



0)+|I A 2 




co+H A 2 



FIG. 1: Density of states pA in the alloy-analogy (solid line) 
and Hubbard-I (dashed line) approximations on the Bethe 
lattice (tA = 1) for various ub and U. 



model with the strong on-site repulsion U. 
irreducible part in this approximation reads 



i 



ui + [l a 



The total 



(54) 



It is well known that the Hubbard-I approximation 
gives an electron spectrum with the opened gap for any 
nonzero repulsion U (Fig. Thus, this approximation 
cannot describe metal-insulator transitions. However, we 
consider it here to estimate its applicability at large U 
from a thermodynamic point of view. 

The interacting single-site Green's function G a (to) is 
calculated using the integration l(T7|) with a noninteract- 
ing density of states p° a {e). The band of particles of a 
given type (a) depends here only on p a (e) and a concen- 
tration of particles of the opposite type, i.e., the transfer 
parameter t s does not have an effect on G a {ui). 

In the large-?/ limit (U ^> t a ), the density of states has 
the simple form 



- U 



(55) 



Both spectral subbands are of the same height for any 
concentration values and are equal to p^.(e) with the 
scaled bandwidth. 



B. Alloy-analogy approximation 

The system with one type of particles frozen on a 
lattice can be mapped onto the problem of the elec- 
tronic structure of simple binary alloys [TtI lil j. The 
alloy-analogy approximation is the single-site solution of 



the binary alloy problem within the coherent potential 
approximation |45j that is exact in infinite dimensions. 
Such a mapping onto the binary alloy is exact for the 
spinlcss Falicov-Kimball model with localized ions. Thus, 
the alloy-analogy approximation exactly describes the 
band formation by the particle hopping (transfer) and 
the interaction with the localized particles. However, in 
the case of the asymmetric Hubbard model this approx- 
imation does not include an effect of the transfer of par- 
ticles of another sort, like the Hubbard-I approximation. 

The self-energy part for the Dyson equation in the 
alloy-analogy approximation reads 



where the single-site Green's function: 
G°M = G a {Lo) 

1 - fla Uif 



(56) 



U) + fl a - J a {u>) LU + [l a - U - Jo-(w) 



(57) 



The Larkin irreducible part S CT 1 = uj + /i CT — S CT is deter- 
mined by the following expression 



1 — n- n 



V + M<r - J a (w) UJ + flu ~ U — J a (uj) 
+ Ja{0j). 



(58) 



This approximation is obtained in our approach when 
the projecting constant (p a and the function R a (uj) in 
lHH} are neglected. 

To compare analytically results of the alloy-analogy 
approximation and the Hubbard-I approximation we can 
examine the equations in the infinite-?/ limit. The self- 
consistency condition gives J a — t\G a on the Bethe lat- 
tice. Thus, in this case, we have the following expression 
for the particle spectrum: 



AA 



1 



H = 7T-T oV4#(l-n*)-(w + ^) 2 



2nt 



(59) 



when \uj + p a \ < 2t a \J\ — n s and zero otherwise. The 
height and width of this spectral band depend on the 
concentration n s and they are scaled proportionally to 
\J\ — jia when n s changes. Therefore, the bandwidth of 
the subbands (for finite U) is wider in the alloy-analogy 
approximation at half-filling Fig. IHA,C) as well as off 
half-filling Fig.[HB,D). In Fig.0 the effect of wider sub- 
bands is illustrated by the dependence of the chemical 
potential p,A on the concentration ha calculated using 
the particle density of states 1)47(1. 



V. ANALYTIC PROPERTIES AT MOTT 
TRANSITION 

We consider here the metal-insulator transition at half- 
filling. It is known that for the Hubbard model at zero 



6 




FIG. 2: The chemical potential as a function of the con- 
centration in the alloy-analogy approximation (solid line) on 
the Bethe lattice compared to the Hubbard-I approximation 
(dashed line) and the atomic limit t a — » (dotted line). 



temperature the transition occurs when Im S CT (0) = 
and the real part has the following low-frequency form 



Re£ CT (w + i0+) = U/2 + (1 - l/Z)u + 0{uj 3 ) (60) 



corresponding to Fermi liquid behaviour, where Z defines 
the quasiparticles weight. Such self-energy gives particle 
densities of states with a quasiparticle peak at the Fermi- 
level that has the same height for different U up to the 
critical value U c [2f| . The gap opens only when Z — > 
and in the metallic phase the relation G a (0)J a (0) = — 1 
is satisfied. 

However, for finite temperatures the Fermi-liquid state 
breaks near the transition point U c . Therefore, S cr (0) 
has a nonzero imaginary part [23] , and the limit Z — > 
cannot indicate the transition. In this case, the density of 
states at Fermi-level continuously tends to zero. At high 
enough temperatures, the Mott transition turns into a 
crossover from a bad metal to a bad insulator ^y] . 

Another limit of the asymmetric Hubbard model (the 
Falicov-Kimball model) at half-filling does not have a 
central quasiparticle peak in the spectrum [23, H3| ■ The 
alloy-analogy, Hubbard-III and GH3 approximations also 
cannot describe the Fermi liquid behaviour at half-filling. 
Hence, these approximations can be applied to the in- 
vestigation of high temperature properties of the system 
when the quasiparticle features disappear. 

Thus, to investigate the continuous transition in the 
asymmetric Hubbard model, we derive some properties 
following from the particle hole-symmetry. 



A. Exact relations between G a (0) , J CT (0) and H CT (0) 
at continuous crossover at half-filling 

In the case of the symmetric noninteracting density of 
states 

p°( £ )=p°(-e) (61) 
leading to the particle-hole symmetry with the properties 
G a {u) = -G a (-w) t JAuj) = -J a {-io) (62) 

of the single-particle Green's function and the coherent 
potential at half-filling, the following relations can be 
proven 



lim 



Mo) 



u^u~ G<j (0) 



-t-oo 



e^(e)de, 



lim (^(OJS-^O) 



1 



(63) 



(64) 



which take place at zero frequency when the continuous 
metal-insulator transition (crossover) is approached from 
below (when the gap opens in the particle spectrum at 
u = 0). 

Proof. Let us consider general properties of the Green's 
functions (G^, J CT ) and the total irreducible part S CT fol- 
lowing from (|62|1 at zero frequency. These functions can 
be written in the Lehmann representation 



— lim 

7T r;^0+ J _ r 



lmF(u' + ir])duj' 



(65) 



This representation shows that due to the symmetry l|62() 
the functions at to — are pure imaginary, and there are 
three possible cases: (i) F(0) = 0; (ii) ReF(0) = 0, 
Imi^O) = const; (iii) the pole of F(uj) in the real axis at 
u = 0. 

An imaginary part of G a (uj ± + ) defines the particle 
spectrum (|47fl . When the gap continuously opens the 
single-particle Green's function tends to zero at u) = 0. 
It means that S~ 1 (0) is pure imaginary 



lim E- 1 {±ir))=±iB a , B a > 0, 



(66) 



and the gap opens (in the non-Fermi liquid state) only 
when Bu continuously increases up to infinity, i.e., the 
self-energy diverges [TJ, |3(j, E3] . 

The coherent potential can be expressed from (|16fl as 
a function of G„ and the total irreducible part S CT : 



Ja 

G„ 



Gl 



(67) 



Using that the noninteracting density of states is sym- 
metric (|61|l and normalized 



p°(e)de 



(68) 
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we have: 



+00 n 



p°(e)de 



2 ' 



and 



GfjJrr — ^„ G a 1 



D 



+ 00 _2 J3 



e p a (e)de 



(69) 



(70) 



from l(T7l) . 

When the critical value U c is approached from below, 
the limit B a — > +00 should be considered. In this limit, 
inserting and (JTOJ into |||)7J) leads to the result (|B3|l . 

To prove (|fi4"|) . we consider three different cases depend- 
ing on behaviour of the function p° (e) at large energies. 

1. If an average kinetic energy per particle is finite, 
the limit 



lim s 2 pl{e) = 



(71) 



is satisfied giving the dependence G a J a ~ B~ 2 (f7T7f> 
at large B a . It means that J CT (0) — > at the tran- 
sition. 

2. When the density of states has the power-law tails 
of order e~ 2 



lim e 2 p®(e) = const > 0, 



(72) 



the limiting behaviour G a J a ~ B" 1 is realized at 
the transition. In this case a finite imaginary part 
of J(t(0) remains when B a — ► 00, and the ratio 
J (T (0)/G (T (0) tends to infinity ^ . 

An example of such a function is the Lorentzian 
density of states describing a lattice with long- 
range hopping along coordinate axes 



n(e 2 +tl) 



(73) 



This density of states always gives a constant value 
for the coherent potential: 



lim J rj (u> ± i?7) = =fit a 



(74) 



3. The last case is when p°(e) ~ £~ n , n £ (1,2) at 
large energy e. The integrals Q68[) and (|7U|) are 
convergent only when n > 1. Such power-law 
tails give the following limiting behaviour G a J a ~ 
l/(B a ) n ~ 1 and J a ~ {B a ) 2 - n at large B a . 



All these cases show that the relation 
lim G CT (0)J ff (0) = - 

[/->£/" 

is satisfied proving the relation (|64|) . 



(75) 



B. Critical value (7 C for the Mott transition within 
GH3 approximation with symmetric density of 
states 



Due to (|63() . we have 

J a (0) = W 2 G a (0)/4 



(76) 



at the Mott transition, where W a is the effective half- 
bandwidth of the unperturbed density of states p® a : 



W a = 2(( + e 2 p° a (e)de 

\J — OG 



1/2 



(77) 



Due to G cr (0)H CT 1 (0) = 1 (|64|l at the transition, insert- 
ing (|?6*j) into (|52"|> yields a set of equations: 



l im U 2 Ga(0) 
u ™- WjG^O) + 2W|G B (0) 



t/ 2 G B (0) 
1*^(0) +2W^G A (0) 



The ratio 



7/ = lim 



Ga(0) 



t^t,- G s (0) 



(78) 



(79) 



(80) 



is excluded from the set of equations, and we have the 
following equation for U c 

U*-(Wl + W 2 )U 2 -3WlW 2 =0 (81) 

with the solution 

Uc _ s l wi + w% + (w\ + wj + uw 2 w 2 y/ 2 _ (82) 

If we put Wb = 0, the expression l|82ll describes the 
Mott-type transition in the alloy-analogy approximation 
(which is exact for the Falicov-Kimball model) giving well 
known results: U c = Wa = 2tA for the Bethe lattice and 
U c = V2tA for the hypercubic lattice. In the case of the 
standard Hubbard model (Wa = Wb = W) , we have the 
result U c — \/~3W corresponding to the Hubbard-III ap- 
proximation at half-filling (see Refs. imilfil for reviews of 
the limits: the standard Hubbard model and the Falicov- 
Kimball model). 

VI. RESULTS AND DISCUSSION 

There are two different pairs of parameters (ua, Pa 
and ns, Pb) m the asymmetric Hubbard model which 
have to be determined from the thermodynamic equi- 
librium. Thus, various thermodynamic regimes may be 
considered when different pairs of the thermodynamic pa- 
rameters are fixed: (pa, Pb), (pa, "s), (n>A, tib), etc. It 
is known that there are possible phase separations at low 
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tem pera tures in the case of the Falicov-Kimball model 
[34l l48l l49| and in the general case of the asymmetric 
Hubbard model Q. Besides the segregated phases, the 
long-range antiferromagnetic-type ordering is possible at 
low temperatures 

Since the investigation of the asymmetric Hubbard 
model is a very complicated problem, we restrict our 
analysis to temperatures higher than the critical temper- 
atures of thermodynamic instabilities. Thus, the problem 
is reduced to the investigation of a band structure and 
metal-insulator transitions at constant particle concen- 
trations. 

The simplest approximations such as the Hubbard-I 
and alloy-analogy approximations describe the band of 
particles generated by the particle hopping and the in- 
teraction with particles of another sort. These approxi- 
mations can describe an effect of the transfer parameter 
tA on the spectrum pb(^>) only via the concentration ua, 
i.e., when the concentrations are fixed the band pb{^>) 
does not depend on tA, and for tA 7^ 0, ts = the ap- 
proximate spectrum pb(w) has a form of two delta-peaks. 
However, the alloy-analogy approximation is exact only 
for itinerant particles in the Falicov-Kimball model when 
another sort of particles is localized. Thus, this approxi- 
mation can give simple reasonable results for dependence 
Pa = Pa(n a , n s ) in the limit of small values of n s or t s . 

The generalization of the Hubbard-III approximation 
(GH3) gives the total irreducible part and the single-site 
Green's function as a functional of the coherent poten- 
tials of both sorts of particles (see equations l]44 p -l]46 p ). 
Unlike the alloy-analogy approximation, GH3 described 
broadening of the band by the interaction with parti- 
cles of another sort, and it gives a finite bandwidth for 
localized particles in the Falicov-Kimball model. In gen- 
eral, the spectral density in the GH3 approximation is 
temperature dependent; such an approximation can be 
applied for systems with different values of the transfer 
parameters (both tA — ts and tA ts) and various 
particle concentrations. However, the approximation has 
some restrictions. In Ref. 0, GH3 was applied to the 
calculation of the chemical potential and spectrum of lo- 
calized particles in the Falicov-Kimball model. It was 
shown that the approximation gives better results when 
the concentration of itinerant particles is low. 

In the symmetric case of half-filling the approximate 
density of states is independent of temperature (Fig.|3J), 
because the projecting coefficients and the coefficients of 
the integral terms in R a {^>) (|46|l are equal to zero. In 
the Falicov-Kimball model at half-filling, the approxima- 
tion (|53|l can be compared with the exact results [2(ll22j . 
We find that the approximation gives a better correspon- 
dence with the exact curves at high temperatures, and 
the best correspondence is for high U when the exact 
results weakly depend on temperature 0. Fig. OJA) 
shows that the density of states of heavy particles B (the 
case of tu < tA) is in the form of a single peak for the 
low interaction U. For the higher interaction strength 
Fig. ElB,C), the spectrum ps(w) has two subbands with 




FIG. 3: Densities of states of particles in the asymmetric 
Hubbard model within the GH3 approximation for various 
values of U at half-filling on the hypercubic lattice. 



peaks which are closer to the centre (u = 0) than in 
the spectrum of light particles pa(u>). This agrees with 
the results obtained for the Falicov-Kimball model when 
is = 0. 

The particle-hole symmetry at half-filling simplifies the 
investigation of the problem. Such symmetry requires 
the divergence (l/u>) of the self-energy at zero frequency 
when the gap opens continuously [30, E3- Thus, the 
limit behaviour l|63|) . (|64|l of the Green's function and the 
coherent potential at the transition is proven analytically, 
which allows us to perform some analytic analysis. 

In Fig. |3 the transition with the opening gap is il- 
lustrated for the asymmetric Hubbard model within the 
GH3 approximation. There arc only continuous transi- 
tions. This agrees with the fact that the Mott transition 
has to be of a continuous type at high temperatures for 
the standard Hubbard model |25l |30|, and in the case 
of the Falicov-Kimball model the gap opens continuously 
with increasing U 0] . The spectra of both types of par- 
ticles have the same bandwidth in the approximation. 
This is because it cannot describe tails of the subbands 
that raise with the temperature increase as it is for lo- 
calized particles in the Falicov-Kimball model [23, l22| . 
Since the effect of temperature broadening of the band is 
not captured, it can be predicted that the approximation 
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FIG. 4: Dependence of the energy band edges on the repul- 
sion strength U at half- filling ua — tib = 1/2 on the Bethe 
lattice: tA = 1, is = 0.5. The solid line is the GH3 approxi- 
mation, the dashed and dotted lines are the alloy-analogy and 
Hubbard-I approximations, respectively. 



tion of the hopping of both types of particles, and the 
critical value U c > U AA . Because of the particle- hole 
symmetry at half-filling the critical value U c can be cal- 
culated analytically in the GH3 approximation and it is 
given by the expression l|82(l . see Fig. [5] The obtained re- 
sult can be applied to various lattices (hypercubic, Bcthc, 
etc.) when the energy is scaled in the following way: 



u: = 



+00 



e 2 P° A (e)de 



-1/2 



(83) 



Uc_ = Uc 

W A 2 

It is common practice to normalize U by a typical ki- 
netic energy (Wa) investigating the Mott transition in 
the standard Hubbard model . In this case, such nor- 
malization shows that the scaled result does not depend 
on a form of the noninteracting density of states. For the 
Lorentzian density of states l|73|) with long-range hopping 
an average kinetic energy per particle tends to infinity 
(Wa — > oo), and as it was noted in Ref. this requires 
an infinite interaction U to drive the system insulating. 



VII. CONCLUSIONS 



1.6 - 

U* ■ 

c 

1.4 - 

1.2 - 

1.0 - 



0.0 0.2 0.4 f /t 0.6 0.8 1.0 

FIG. 5: Scaled critical value U* = U C /W A EBt as a function 
of t b ft a within the GH3 approximation at half-filling. 

underestimate the critical value U c . 

The energy band edges as a function of U within vari- 
ous approximations on the Bethe lattice are displayed in 
Fig. 2| In the Hubbard-I approximation the spectrum is 
always split having the gap and two subbands. The alloy- 
analogy approximation is exact for the band of itinerant 
particles when ts = and it gives the critical interaction 
constant U AA — 2t A - The GH3 approximation describes 
the band which is broadened by simultaneous contribu- 



The generalization of the Hubbard-III approximation 
is obtained for the asymmetric Hubbard model using the 
equation of motion approach. This method combines in 
a unified framework the description of the band structure 
of the Falicov-Kimball model and the standard Hubbard 
model. In general, the approach can be used to calculate 
the particle spectrum in the system with the different 
particle concentrations at various temperatures, and it 
gives nontrivial results for the spectrum of localized par- 
ticles in the Falicov-Kimball model. 

The self-consistency condition connecting mutually the 
single-site Green's function and the coherent potential 
gives the simple result (I63fl when the metal-insulator 
transition occurs at half-filling on the lattice with a sym- 
metric noninteracting density of states. As a result, we 
have the expression (|82|) in the GH3 approximation for 
the critical value of U that provides a universal solution 
for various lattices and different values of the transfer 
parameters. 

The approximation has some restrictions. Thus, to de- 
scribe the system at half-filling with the proper tempera- 
ture dependence or the quasi particle peak, the approach 
needs further improvement in calculating the mass oper- 
ator by including higher order corrections. 
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